Measurement of the Casimir force between a gold sphere and sihcon surface with 

nanoscale trench arrays 
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We report measurements of the Casimir force between a gold sphere and a silicon surface with 
an array of nanoscale, rectangular corrugations using a micromechanical torsional oscillator. At 
distance between 150 nm and 500 nm, the measured force shows significant deviations from the 
pairwise additive formulism, demonstrating the strong dependence of the Casimir force on the 
shape of the interacting bodies. The observed deviation, however, is smaller than the calculated 
values for perfectly conducting surfaces, possibly due to the interplay between finite conductivity 
and geometry effects. 
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The Casimir force is the interaction between neutral 
conductors that can be understood as resulting from the 
alteration of the zero point energy of the electromagnetic 
field in the presence of boundaries [ij. For two perfect 
metallic planar surfaces, the force is attractive and is 
given by Fc = 7r^?icv4/240z*, where c is the speed of light, 
fi is the Planck's constant/27r, z is the separation between 
the plates and A is the area of the plates. There exists 
a close connection between the Casimir force between 
conductors and the van der Waals (vdW) force between 
molecules. For the former, the quantum fluctuations are 
often associated with the vacuum electromagnetic field, 
while the latter commonly refers to the interaction be- 
tween fluctuating dipoles. In simple geometries such as 
two parallel planes, the Casimir force can be interpreted 
as an extension of the vdW force in the retarded limit. 
The interaction between molecules in the two plates is 
summed to yield the total force. However, such summa- 
tion of the vdW force is not always valid for extended 
bodies because the vdW force is not pairwise additive. 
The interaction between two molecules is affected by the 
presence of a third molecule. One important characteris- 
tic of the Casimir force is its strong dependence on geom- 
etry [31. The Casimir energy for a conducting spherical 
shell [3| or a rectangular box [1, has been calculated to 
have opposite sign to parallel plates. Whether such ge- 
ometries exhibit repulsive Casimir forces remains a topic 
of current interest |6|. 

In recent years, there have been a number of precision 

measurements of the Casimir force 0, i, i, Efl, in m, 

13, 14, These experiments yield agreement with the 
theoretical calculations to accuracies of better than one 
percent when nonideal behavior of the metallic surfaces 
[la . [iTl [l8 | are taken into account. The vast majority 
of force measurements were performed between a sphere 
and a flat plate, two flat plates or two cylinders. For these 
simple geometries, the Casimir force is not expected to 
show significant deviations from the pairwise additive ap- 




FIG. 1; (a) Cross section of rectangular trenches in silicon, 
with periodicity of 400 nm and depth of 0.98 /xm (sample B). 
(b) Top view of the structure, (c) Schematic of the exper- 
imental setup (not to scale) including the micromechanical 
torsional oscillator, gold spheres and silicon trench array, (d) 
Measurement scheme with electrical connections. Excitation 
voltages Vaci and Vac2 are applied to the bottom electrodes. 



proximation (PAA) at small separations. There has only 
been one experiment that involved surfaces of other ge- 
ometries, where the Casimir force is measured between a 
sphere and a plate with small sinusoidal corrugations [l^ . 
While this measurement shows deviations from PAA, the 
interpretation of the deviation is still controversial. For 
example, the Casimir force for sinusoidal corrugations on 
perfect conductors was calculated without the assump- 
tion of PAA [i^l . The deviation from PAA is found to 
be strong only when the ratio of the separation to the 
periodicity of corrugation A is large. In ref. [31, A is not 
small enough for deviations from PAA to be significant. 
It is suggested that lateral movement of the two surfaces 
may be able to account for the deviations 0| . 

In this letter, we report measurements of the Casimir 
force between nanostructured silicon surfaces and a gold 
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sphere. One of the interacting objects consists of a sil- 
icon surface with nanoscale, high aspect ratio rectangu- 
lar corrugations. The other surface is a gold-coated glass 
sphere attached onto a micromechanical torsional oscilla- 
tor. Lateral movements of the surfaces are avoided by po- 
sitioning the corrugations perpendicular to the torsional 
axis. The Casimir force gradient is measured from the 
shifts in the resonant frequency of the oscillator at dis- 
tance between 150 nm and 500 nm. Deviations of up to 
20% from PAA is observed, demonstrating the strong ge- 
ometry dependence of the Casimir force. The measured 
deviation is, however, about a factor of 2 smaller than de- 
viations expected for perfectly conducting surfaces [22 1. 

Figure [1^ shows a cross section of an array of rect- 
angular corrugations with period of 400 nm (sample B) 
fabricated on a highly p-doped silicon substrate. Two 
other samples, one with period 1 ^m (sample A) and the 
other with a flat surface, are also fabricated. The fab- 
rication procedure started with a layer of silicon oxide 
(0.2 fj,m) deposited onto a blank silicon wafer by chemi- 
cal vapor deposition. Lithography was performed with a 
deep ultra-violet stepper followed by reactive ion etching 
to transfer the pattern from the photoresist to the silicon 
oxide. Trenches with depths t — 2a {^1 /im) were then 
created using deep reactive ion etching using the oxide 
as etch mask. A continuous etch and deposit recipe was 
used to yield smooth side walls at 90.3 ° and 91.0 ° to the 
top surface respectively for samples A and B. Residual 
hydrocarbons were removed using an oxygen plasma etch. 
Finally, the oxide etch mask is removed using hydroflu- 
oric acid. In order to ensure that the optical properties 
of the silicon are identical, all samples were fabricated 
on the same wafer and later diced into 0.7mm by 0.7mm 
pieces for the force measurement. 

The geometry of nanoscale, rectangular trenches was 
chosen because the Casimir force on such structures are 
expected to exhibit large deviations from PAA. We con- 
sider the interaction between the trench array and a par- 
allel flat surface at distance z from the top surface of the 
trenches. In the pairwise additive picture, this interac- 
tion is a sum of two contributions: the volume from the 
top surface to the bottom of the trench and the volume 
below the bottom of the trench. The latter component 
is negligible because the distance to the other surface is 
more than 1 ^m, larger than the distance range at which 
Casimir forces can be detected in our experiment. For 
a trench array of 50% duty cycle, the former component 
yields exactly half of the interaction between two flat 
surfaces Ffiat regardless of the periodicity because half 
of the material is removed In practice, the trench ar- 
rays are created with duty cycle close to but not exactly 
at 50%. Under PAA, the total force is equal to pFfiat, 
where p is the fraction of solid volume. The calculation of 
the Casimir force in such corrugated surfaces, in contrast, 
is highly non-trivial. While perturbative treatments [l^ 
are valid for smooth profiles with small local curvature. 



they are impractical for the deep, rectangular corruga- 
tions. Using a different approach based on path inte- 
grals, Biischer and Emig 22| calculated the Casimir force 
for the corrugated geometry made of perfect conductors. 
Strong deviations from PAA were obtained when the ra- 
tio z/X is large. In the limit when A goes to zero, the 
force on a trench array approaches the value between flat 
surfaces, leading to deviations from PAA by a factor of 
2. Such large deviations occur because the Casimir force 
is associated with confined electromagnetic modes with 
wavelength comparable to the separation between the in- 
teracting objects. When A << z, these modes fail to 
penetrate into the trenches, rendering the Casimir force 
on the corrugated surface equal to a flat one. 

We measure the gradient of the Casimir force on the 
silicon trench arrays using a gold-coated sphere attached 
to a micromechanical torsional oscillator [ld|. The oscil- 
lator consists of a 3.5^m thick, 500/im square silicon plate 
suspended by two torsional rods. As shown in Fig{T];, two 
glass spheres, each with radius R of 50 ^m, are stacked 
and attached by conductive epoxy onto the oscillator [l3| 
at a distance of 6 = 210 fim from the rotation axis. The 
large distance 200 /xm) between the oscillator plate 
and the corrugated surface ensures that the attraction be- 
tween them is negligible and only the interaction between 
the top sphere and the corrugated surface is measured. 
Before attachment, a layer of gold with thickness 4000 A 
is sputtered onto the spheres. Two electrodes are located 
between the plate and the substrate. Torsional oscilla- 
tions in the plate are excited when the voltage on one of 
the electrodes is modulated at the resonant frequency of 
the oscillator (/o = 1783 Hz, quality factor Q = 32000). 
For detecting the oscillations, additional ac voltages at 
amplitude of 100 mV and frequency of 102 kHz is applied 
to measure the capacitance change between the top plate 
and the electrodes. A phase locked loop is used to track 
the shifts in the resonance frequency [l0| as the sphere 
approaches the other silicon plate through extension of a 
closed- loop piezoelectric actuator. As shown in Fig.[T};, 
the movable plate is positioned so that its torsional axis 
is perpendicular to the trench arrays in the other sili- 
con surface. Such an arrangement eliminates motion of 
the movable plate in response to lateral Casimir forces 
21 , 2^ because the spring constant for translation along 
the torsional axis is orders of magnitude larger than the 
orthogonal direction in the plane of the substrate. 

To prepare the silicon surface [l^ for force measure- 
ment, hydrofluoric acid is used to remove the native oxide 
on the surface of the silicon chip. The hydrofluoric acid 
also passivates the silicon surface to temporarily prevent 
oxide formation at ambient pressure. In the next step, 
the silicon chip is baked at 120 °C for 15 minutes to elim- 
inate residual water that might have accumulated in the 
trenches. The silicon chip is then positioned to within a 
few micrometers from the gold sphere and the chamber 
is immediately evacuated to a pressure of 10~^ torr by a 
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FIG. 2: Gradient of the electrostatic force {V = Vb + 300mV) 
on the fiat silicon surface (solid circles) and corrugated sili- 
con structure (hollow squares). The solid line is a fit using 
Eq. [2] for the flat surface. The dashed line is a fit using the 
force gradient from finite element analysis for the corrugated 
structure. Inset: The space between the corrugated structure 
and a flat surface is divided into triangular mesh to solve the 
Poisson equation in 2D (2 = 200 nm) . 



dry roughing pump and a turbo pump. 

For small oscillations where nonlinear effects can be 
neglected, the resonant frequency shift A/ of the oscil- 
lator is proportional to the gradient of the force F'{z) 
between the surfaces: 



A/ = CF'(z), 



(1) 



where C = — 6^/87r^/o/ and / is the moment of inertia of 



the oscillator. The distance z is given by zq ■ 



^piezo 



where zq is the initial separation between the two sur- 
faces, Zpiezo is the extension of the piezoelectric actuator 
and b9 is the modification of the separation due to ro- 
tation of the top plate to angle 9. To calibrate zq and 
the proportionality constant C, a dc voltage V is ap- 
plied to the silicon plate with the gold sphere electrically 
grounded. The electrostatic force between a sphere and 
a flat plate is given by: 



2neo{V - Vo) Mh(a) - n coth(na) 



sinh(na) 



(2) 



where cosh a = 1-1- d/R and Eq is the permittivity of 
vacuum. The residual voltage Vq is measured to be ^ 
—0A3V by determining V at which A/ attains minimum 
at fixed z. For 100 nm < z < 2 iim, Vq is found to vary 
by less than 3 mV. In Fig. [21 the solid line represents a 
fit to the measured electrostatic force gradient for a fiat 
surface. The contribution of the Casimir force gradient 
to the measured frequency shift (< 4% at the smallest 
z) has been subtracted. By averaging the fitted values 
for six sets of data with V ranging from Vq + 245 mV to 
Vo + 300 mV, C is determined to be 628 ±5mN-^ s-^. 
The voltages are chosen to be larger than Vq to avoid 
depleting the surface of the p-doped silicon with charge 
carriers. A similar calibration procedure is performed on 



the silicon trench arrays (dashed line in Fig. [2]). Since 
there is no analytic expression for the trench geometry, 
the electrostatic force is calculated by solving Poisson's 
equation in 2D using finite element analysis. As shown 
in the inset to Fig. [21 the boundary conditions are set by 
maintaining a fixed potential between the trench array 
and a fiat surface, with the volume between them divided 
into N > 10, 000 triangles. Since R >> z, the proximity 
force approximation Fsc — 2'nREfc can be used to obtain 
the force F^c between a sphere and a corrugated surface 
from the electrostatic energy Efc between a fiat surface 
and a corrugated surface. To ensure that N is sufficiently 
large, we checked that the calculated force varies by less 
than 0.1% even when N is doubled. 

The dimensions of the trench array used in the cal- 
culation were obtained from the scanning electron mi- 
crographs: the fraction of solid volume p is determined 
from the percentage of bright pixels in the top view sim- 
ilar to Fig.[T]3. Ten pictures at different locations with 
area 30 /im by 50 /xm for sample A and 20 /xm by 35 /xm 
for sample B were used to calculate pa and pb to be 
0.478 ± 0.002 and 0.510 ± 0.001 respectively. The cross 
sectional view yields the trench depth [tA = 0.98 fim and 
Ib = 1.07 /im). As shown in Fig. la, while we have op- 
timized the fabrication process so that the top corner of 
the trenches has a sharp rectangular shape, the bottom 
sections show certain degree of rounding. In all analysis 
described here, the trenches are assumed to have perfect 
rectangular shape. The validity of such approximation 
is justified by the insensitivity of the calculated electro- 
static force on the depth of the trenches. Varying the 
depth of the trenches by 10% produces less than 0.01% 
change in the calculated force. 

By setting V equal to Vq, the Casimir force gradient 
F^ /;at is measured between the gold sphere and a flat sil- 
icon surface (solid circles in Fig. [3^) obtained from same 
wafer on which the corrugated samples A and B were 
fabricated. The main source of uncertainty in the mea- 
surement (~ 0.64pN/im^^ at z = 300nm) originates from 
the thcrmomcchanical fluctuations of the micromechan- 
ical oscillator. As the distance decreases, the oscillation 
amplitude is reduced to prevent the oscillator from enter- 
ing the nonlinear regime [13] ■ At distances below 150 nm, 
the oscillation amplitude becomes too small for reliable 
operation of the phase-locked loop. In Fig. [3^, the line 
represents the theoretical force gradient between the gold 
sphere and the flat silicon surface, including both the fi- 
nite conductivity and roughness corrections. Lifshiftz's 
expressions [l3, [3, [l3] are used to take into account the 
finite conductivity. For the gold surface, tabulated values 
of the optical properties [2J] were used. For the silicon 
surface, the tabulated values were further modified by the 
concentration of carriers (2 x lO^^cm"'^) determined from 



the dc conductivity of the wafer (0.028 cm) [IJ]. Using 
an atomic force microscope (NTMDT), the main contri- 
bution to the roughness is found to originate from the 
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FIG. 3: Measured Casimir force gradient between the same 
gold sphere and (a) a flat sihcon surface, F^ ji^^, (b) sample A, 
F'c,A = l^m) and (c) sample B, F^ g (A = 400nm). In (a), 
the line represents the theoretical Casimir force gradient in- 
cluding finite conductivity and surface roughness corrections. 
In (b) and (c), the lines represent the force gradients expected 
from PAA {pF'^ji^i). (d) Ratio p of the measured Casimir 
force gradient to the force gradient expected from PAA, for 
samples A (A/a = 1.87, hollow circles) and B (A/a = 0.82, 
solid squares) respectively. Theoretical values [23] for per- 
fectly conducting surfaces are plotted as the solid (A/a = 2) 
and dashed lines (A /a = 1). 



gold surface 4 nm rms) rather than the silicon wafer 
{'^ 0.6 nm rms) 
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The roughness correction 
taken into account by the geometrical averaging method 

The Casimir force gradients F'^ ^ and F'^ g between the 
same gold sphere and the corrugated samples A and B 
were then measured and plotted in Figs. [3)3 and [21; . As 
described earlier, under PAA, the forces on the trench 
arrays (where z is measured from the top of the corru- 
gated surface) are equal to the force on a flat surface 
multiplied by the fractional volumes pA and pb- The 
solid lines in Figs. [5Jd and[5j3 represent the correspond- 
ing force gradients, pA F^ ^^^^ and pb F^ ^^^^ respectively. 
Measurement of the force gradient was repeated three 
times for each sample, yielding results that are consistent 
within the measurement uncertainty. To analyze the de- 
viations from PAA, the ratios pA — F'^ aIva F'^ jj^j and 
Pb = F'^ bIpb Pc flat ^^'^ plotted in Fig. [5Ji. The ratio p 
equals one if PAA is valid. For sample A with A/a = 1 .87, 
where a is half the depth of the trenches, the measured 
force deviates from PAA by ~ 10%. In sample B with 
X/a ~ 0.82, the deviation increases to ^ 20%. For 150nm 
< z < 250 nm, the measured Casimir force gradient in 
both samples show clear deviations from PAA. At larger 
distances, the uncertainty increases considerably as the 
force gradient decreases. 

We compare our experimental results on silicon struc- 



tures to calculations by Biischer and Emig [24I on perfect 
conductors (solid and dashed lines in Fig. [3li). In this 
calculation, the Casimir force between a flat surface and 
a corrugated structure with p = 0.5 was determined for 
a range of A/a using a path integral approach. Since 
R >> z, the proximity force approximation allows a di- 
rect comparison of our measured force gradient using a 
sphere and the predicted force that involved a flat sur- 
face. The measured deviation in sample B is larger than 
sample A, in agreement with the notion that geometry 
effects become stronger as A/a decreases. However, the 
measured deviations from PAA are smaller than the pre- 
dicted values by about 50%, significantly exceeding the 
measurement uncertainty for 150nm < z < 250nm. Such 
discrepancy is, to a certain extent, expected as a result 
of the interplay between finite conductivity and geometry 
effects. The relatively large value of the skin depth in sil- 
icon (~ 1 1 nm at wavelength of 300 nm) could reduce the 
deviations from PAA. So far, exact computation of the 
Casimir force including both finite conductivity and ge- 
ometry effects for strongly deformed rectangular trenches 
has not yet been performed. It is a nontrivial problem 
that warrants further theoretical analysis. 

In conclusion, we measured the Casimir force gradient 
between a spherical gold surface and a silicon surface 
with an array of deep, rectangular trenches. Deviation 
from the PAA is clearly observed and found to increase 
as A/a decreases, consistent with the strong boundary 
dependence of the Casimir force. However, the measured 
deviation is smaller than the calculated value on perfect 
conductors. 
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